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ABSTRACT 

We present the analysis of archival X-ray observations of the supernova remnant (SNR) G21.5— 0.9. 
Based on its morphology and spectral properties, G21.5— 0.9 has been classified as a Crab-like SNR. For 
that reason, it was chosen as a Chandra calibration target. In their early analysis of part of these calibra- 
tion data, Slane et al. (2000) discovered a low-surface-brightness, extended emission. They interpreted 
this component as the blast wave formed in the supernova (SN) explosion. XMM-JVewton observations 
(Warwick et al. 2000) revealed the non-thermal nature of this emission suggesting that it is instead an 
extension of the synchrotron nebula. In this paper, we revisit the Chandra analysis using new calibration 
data, improving the statistics by a factor of 2. We also include ROSAT and ASCA observations. Our 
analysis confirms the non-thermal nature of the extended emission. ACIS-S images indicate that this 
component is not limb-brightened, and it shows knotty structures and a bright filament 2' north of the 
center. We find no evidence of line emission from any part of the remnant. We can reject a collisional 
equilibrium ionization (CEI) thermal model at solar abundances, and non-equilibrium ionization (NEI) 
models (such as a plane-parallel shock model with different ionization ages and constant-temperature, or 
an NEI model with a single ionization age and a constant temperature) . The entire remnant is best fitted 
with a power law model with a photon index steepening away from the center. The total unabsorbed 
flux Fx (0.5 - 10 keV) is 1.1 x 10~^° erg cm~^ s~^ with an 85% contribution from the 40" radius core. 
Timing analysis of the High-Resolution Camera (HRC) data failed to detect any pulsations. We put a 
16% upper limit on the pulsed fraction. We derive the physical parameters of the putative pulsar and 
compare them with those of other plerions (such as the Crab and 3C 58). 

Subject headings: ISM: individual (G21.5— 0.9) - stars: neutron - supernova remnants - X-rays : ISM 



1. INTRODUCTION 

Crab-like supernova remnants have a center-brightened 
nebula, referred to as a plerion, often seen in the radio 
and X-ray wavelengths, and believed to be powered by an 
embedded pulsar. Typified by the Crab nebula, they have 
non-thermal spectra and lack the evidence of a SNR shell 
at all wavelengths. In the radio, they are highly polarized 
and have a flat power law spectral index {a '^ 0—0.3; Si, 
^ z/~"). In X-rays, they have a hard photon index (F ^ 
2; where F = a+l). Out of the ^ 220 Galactic supernova 
remnants (SNRs), only ~ 10 candidates are classified as 
Crab-like supernova remnants (Green 2000). 

G21.5— 0.9 is a plerionic SNR. In the radio, it has 
been imaged with the Very Large Array (VLA) at 5 GHz 
(Becker and Szymkowiak, 1981), with the NRO Millime- 
ter Array at 22.3 GHz, and with the Effelsberg 100-m tele- 
scope at 32 GHz (Morsi & Reich 1987). These observations 
show a centrally peaked highly polarized emission, of ra- 
dius ^ 40", having a flat radio spectral index (a ~ 0) and a 
flux density of 6 Jy at 1 GHz. The 22.3 GHz map obtained 



with the NRO (Fiirst et al. 1988) very closely resembles 
the 5 GHz map of Becker and Szymkowiak (1981, Fig. 1), 
and exhibits an ellipse of 1'.5 major axis, an axial ratio of 
0.8 and a position angle ~35°. 21 cm absorption data ob- 
tained with the VLA place the SNR at a distance between 
4.8 and 9.3 kpc (Becker, & Szymkowiak 1981; Davelaar, 
Smith, & Becker 1986). Assuming a distance of 5 kpc, the 
linear size of the radio plerion is 2.2 pc x 1.3 pc. Fiirst 
et al. (1988) found, in addition to the diffuse centrally 
peaked component, axisymmetric filaments which they in- 
terpret as signatures of jets from a central pulsar. Infrared 
observations of G21.5— 0.9 with ESA's Infrared Space Ob- 
servatory (ISOCAM; Gallant & Tuffs 1998) show that the 
plerion has a morphology similar to the 22 GHz map of 
Fiirst et al. (1988). 

X-ray observations of G21.5— 0.9 point unambigu- 
ously to a centrally peaked non-thermal X-ray plerion. 
G21.5— 0.9 was observed with the non-imaging medium- 
energy (ME) instrument on-board EXOSAT. Analysis of 
these data indicates a power law spectrum with a pho- 
ton index of F = 1.72 ± 0.12, somewhat flatter than 
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that of the Crab (Davelaar et al. 1986). A subsequent 
short observation (2816 s) by GINGA in the 2-10 keV 
range could not distinguish between a power law model 
with a photon index of F = 1.88 ± 0.05 and a ther- 
mal bremsstrahlung model with a very high temperature 
kT = 15.4 ± 1.4 keV (Asaoka & Koyama 1990). While 
the thermal bremsstrahlung yields an acceptable fit to the 
data, not only is the temperature too high but the upper 
limit on the equivalent width of the Fe-K line emission 
(less than 0.17 keV) is too low to accommodate the ex- 
pected value of emission from a plasma at cosmic abun- 
dance in ionization equilibrium. Although the proper- 
ties of G21.5— 0.9 suggest the presence of a pulsar, X-ray 
and radio imaging show no indication of a point source 
(Becker & Szymkowiak 1981, Kaspi et al. 1996), and all 
the searches for pulsation at any wavelength have proved 
inconclusive. G21.5— 0.9 was recently observed with the 
Chandra X-ray observatory as a calibration target. Slane 
et al. (2000) discovered a low-surface-brightness compo- 
nent extending to a radius of ~2' from the central core, 
using a 34 ksec exposure with ACIS-S and a 30 ksec ex- 
posure with HRC-1. The authors interpret this newly dis- 
covered component as the first evidence for the expanding 
ejecta and blast wave formed in the initial explosion. 

In this study we have carried out a combined analysis of 
Chandra, ROSAT, and ASCA datasets. We focus on the 
Chandra data which provide the best angular resolution 
to date. We compile all publicly available observations 
by Chandra and present here the imaging, detailed spec- 
tral and timing analysis of a net 72 ksec exposure with 
ACIS-S and 76 ksec observation with the HRC. This al- 
lowed us to examine in more detail the spectral variations 
across the remnant and put a more stringent upper limit on 
the pulsed fraction from the yet undetected pulsar. With 
the longer exposure with ACIS, we detect the low-surface- 
brightness, extended component (referred to as the 'shell' 
in Slane et al. 2000) out to a radius of 2'. 5, and rule out its 
thermal interpretation. We conclude that the spectrum of 
the entire remnant is dominated by non-thermal emission. 
The broadband images, the hardness ratio map, and the 
spectral properties of the remnant indicate that the entire 
X-ray emission could be due to synchrotron radiation from 
high-energy electrons injected by the 'hidden' pulsar in the 
form of either a wind or jets interacting with the surround- 
ing medium. We discuss our interpretation in comparison 
with the high-resolution 22.3 GHz observations and com- 
pare the parameters of G21.5— 0.9 plerion to the well stud- 
ied Crab and 3C 58 plerions. While G21.5— 0.9 bears sim- 
ilarities to plerions, it has unique properties that makes it 
an intriguing object. To date, it is the only plerion whose 
X-ray size is bigger than in the radio. Its study is im- 
portant to shed light on plerions of the 'second kind', a 
subclass of plerions having a spectral break at millimeter 
wavelengths (Woltjer et al. 1997). The paper is organized 
as follows: In the following section, we describe the obser- 
vations. We describe our imaging analysis in § 3, and the 
corresponding spectral analysis in § 4. In § 5, we present 
the results of our timing analysis. Finally, we discuss our 
results in the last section. 



2. OBSERVATIONS 

2.1. Chandra 

Because of its power law spectrum, small size and 
relative brightness (its unabsorbed flux is ~4.3xl0~^^ 
ergs cm~^ s'-^ between 0.4 and 2.0 keV), G21.5— 0.9 was 
selected as one of the calibration targets of the Chandra 
observatory (Weisskopf, O'dell, & van Speybroeck 1996). 
As such, the SNR was observed in many configurations 
chosen by the instrument teams. In the following anal- 
ysis, we have made use of these different configurations 
to probe different aspects of the remnant. Because of the 
CTI problems^ with the front-illuminated chips, we have 
not used any data from the ACIS-I, even for imaging pur- 
poses (the spatial resolution of about 1" of the damaged 
CCDs is still intact) because the CTI problems make it 
hard to perform accurate energy cuts on the data. 

For the spectral analysis, we used six observations (Ob- 
sID 159, 1230, 1433, 1716, 1717, and 1718) correspond- 
ing to a total exposure time of 71.6 ks and obtained with 
ACIS-S3 (back illuminated chip) for which the most accu- 
rate calibration exists at this point. High- resolution imag- 
ing and timing analysis were done using both the HRC-I 
and HRC-S and their unparalleled angular resolution. 

2.2. ROSAT and ASCA 

G21.5— 0.9 was observed with the ROSAT observatory 
using both the Position Sensitive Proportional Counter 
(PSPC) and the High-Resolution Imager (HRI). Both ob- 
servations are public^, and although they are quite short 
(7.2 ks for the PSPC and 4.7 ks for the HRI), we use the 
PSPC data in the following analysis to better constrain the 
value of the column density and the HRI data to compare 
with the images extracted from the Chandra datasets. 
G21.5— 0.9 was also observed by both the Solid-State 
Imaging Spectrometer (SIS) and the Gas Imaging Spec- 
trometer (CIS) on-board ASCA and the data are also 
available through the HEASARC database. The exposure 
time is 28 ksec with the CIS and 24 ksec with the SIS. Be- 
cause of the relatively poor angular resolution (3') of the 
ASCA X-ray telescopes, none of the fine spatial structures 
of G2 1.5— 0.9 can be studied with neither the CIS nor the 
SIS. ASCA advantages reside in the large number of counts 
collected and the GIS timing capabilities (a timing reso- 
lution of 6.1xl0~^s for the high-bit rate mode in which 
the data were taken). We have used the spectral analy- 
sis of the ASCA datasets as a measure of the "average" 
spectrum from the remnant to complement that extracted 
from the ROSAT PSPC, and in combination with that of 
Chandra, which is the focus of this paper. 

3. IMAGING 

We have screened the ACIS-S data by running /select on 
the primary events files. Only events with standard ASCA 
grades 02346 are retained. We rejected columns which lie 
on the edges of read-out node areas and where event se- 
lection will not operate correctly since the charge cloud 
is split between 2 different read-out nodes. We rejected 
the hot columns and pixels that arc listed on the Chan- 
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dra X-ray Center web site^. We have extracted images of 
G21.5— 0.9 in different energy bands. The images are cen- 
tered at RA (J2000) = 18''33'"33".50, and DEC (J2000) = 
-10°34'6".6. In Fig. 1, we show the soft (0.1-6 keV) and 
hard (6-10 keV) band images, smoothed with a Gaussian 
with (7 = 3". While the extended component is evident in 
the soft band image, only the inner core is left in the hard 
band. The inner core has an angular radius of 40", similar 
to the radio plerion. We refer to this component as the in- 
ner core throughout the paper. The extended component 
extends out to a radius of 150", shghtly bigger than pre- 
viously reported (120", Slane et al. 2000), and confirmed 
by XMM-Newton (Warwick et al. 2000). We refer to 
this component (50"-150") as the low-surface-brightness 
extended component throughout the paper. 

In Fig. 2, we show the hardness ratio map using the 
ratio: (2.4-10) keV over (0.5-2.4 keV). The image is 
smoothed with a Gaussian with a = 3". It is evident 
from the image that the central core is harder than the 
extended component. Moreover, the extended component 
is harder in the northern quadrant, with brighter knots 
and filamentary structures extending from the core out to 
the outer part of the SNR. This result is consistent with 
the detailed spectral analysis described below (§ 4.1). 

In Fig. 3, we show the ACIS-S3 image in the 0.5-2.4 
keV band with the ROSAT HRI contours overlayed. We 
smooth the Chandra image with a Gaussian with a ~ 6" , 
to match the angular resolution of the HRI. The HRI im- 
age reveals the ~40" radius inner core detected in the ra- 
dio. The extended component is detected with ACTS out 
to a radius of 150". The point-like source south-west of 
the plerion coincides with an emission line star, SS 397. 

G21.5— 0.9 was observed five times with the HRC-I for 
a total of more than 100 ks. Unfortunately, two of those 
observations were unusable because of processing or recon- 
struction problems. We have used the three remaining ob- 
servations and produced an image using a mosaic program 
which takes into account the difference in pointing. The 
background is subtracted using a script, screen^hrc, made 
available by the HRC Calibration team^. This procedure 
flags events with a non- X-ray charge cloud distribution and 
removes about 30%-40% of the background events while 
removing only a few percent of real events. We use this 
script and find that 33% of events are cut, leaving a total of 
about 3 X 10^ events. We do not correct for exposure. The 
result is shown in Fig. 4. The image has been convolved 
with an elliptical top hat filter whose size varies to have a 
minimum of 16 counts under the filter. The pixel size is 
0".15. We find that the center of the HRC emission is lo- 
cated at RA (J2000) = 18'^33'"33".51 (±0.01), Dec (J2000) 
= -10°34'8".47 (±0.15). We did not measure the spectral 
difference across the different parts of the HRC-I emission 
because of the lack of spectral response in the HRC. 

In Fig. 5, we show the background-subtracted radial pro- 
file of the plerion obtained with HRC and ACIS. We choose 
the observations with minimal offsets from the telescope 
axis, for which a 0".5-radius aperture encircles 90% of the 
point-source counts. We note that while the ACIS profile 
is shown out to a radius of ~ 100", the HRC profile ex- 
tends to a radius of 40" only due to the high background 
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at larger radii. In the inset, we show the HRC profile with 
0".25 bins along with the expected point-spread function 
(dashed line). This shows that the brightest inner core is 
not a point source and that the brightness drops sharply 
beyond a radius of ~ 0".5. We use this profile later (§ 6.2) 
to infer the luminosity of the putative pulsar in G21.5— 0.9. 

4. SPECTRAL ANALYSIS 
4.1. ACIS-S 

In the following, we present the spectral analysis of the 
SNR and examine the spectral variations on the arcsec- 
ond scale. For the inner 40" core, we extract spectra from 
concentric rings of width 5" and centered at RA (J2000) 
= 18h33™33".5, Dec (J2000) = -10°34'6"6. We also ex- 
tract the spectrum of the brightest 1" radius central part 
of the image. For the extended component, we extract a 
spectrum of a ring extending from 50" to 150" — the latter 
defines the outer boundary detected with Chandra. 

The background is extracted from the same detector re- 
gion on the S3 chip using an observation of G21.5— 0.9 ob- 
tained with the four ACIS-I chips and two ACIS-S chips 
activated. This method allows one to account for both 
the chip and sky background. The background-subtracted 
count rates in the 0.5-10 keV energy range are 2.544 ± 
0.008 ct s-i and 0.506 ± 0.005 ct s'^ from the inner 
core and the extended component, respectively. We pro- 
cessed the data as described in the imaging section. We 
used CIAO vl.l mkarf tool in order to generate the corre- 
sponding ancillary response files. The innermost 1 "-radius 
core has a very hard spectrum and is best described by a 
power law model with F = 1.4 ± 0.1 {Nh = 2.1 ± 0.1 x 
10^^ cm~^), and an unabsorbed fiux in the 0.5-10 keV of 
3.5 X lO"""^^ ergs cm~^ s~^. We subsequently fit all the 
spectra simultaneously with a power law and find that the 
power law index steepens with increasing radius. When al- 
lowing Nh to vary, we get an acceptable fit (x^ — 0.68, i^ = 
4451) with the photon index increasing from 1.62 ± 0.05 at 
the inner 5" circle {Nh = 2.37 ± 0.08 x 10^^ cm^^) ^ 2.36 
± 0.07 at the outer ring {Nh = 1.83 ± 0.07 x 10^^ cm^^). 
The variation of Nh across the plerion could be an artifact 
of the spectral fitting. When tying Nh, we get acceptable 
fits for all the spectra {xl = 0.71; ly = 4459), with Nh= 
2.24 ± 0.03 X 10^^ cm~^ (a value consistent with that 
found by Warwick et al. 2000), and F varying from 1.53 
± 0.03 (inner 5" core) to 2.73 ± 0.04 (the extended com- 
ponent). All errors throughout the paper are at the 90% 
confidence level. We shall accept this value as the Nh for 
the SNR, with an uncertainty we shall determine by the 
range of values we find for fits to other datasets below. 
In Fig. 6, we plot the corresponding photon index as a 
function of radius (with Nh = 2.24 x 10^^ cm~^). The 
steepening of the photon index with increasing radius has 
been noted by Slane et al. (2000) when fitting the inner 40" 
core. Slane et al. (2000) did not however rule out a ther- 
mal model for the extended component, and interpret that 
component as the SNR shell or shocked ejecta. With this 
longer Chandra observation, we find that thermal models 
are unlikely, as discussed further below. 

In Table 1, we summarize the result of fitting the ex- 
tended component with thermal and non-thermal mod- 



els, in the 0.6-8 keV energy range. A power law 
model gives the best fit; with the following parame- 
ters: Nh = l-^it^Qm xlO^^ cm-2, T = 2.36to;o6, and 
X^ = 1.0 (i^=658, see Table 1). When freezing Nh 
to 2.24 X 10^^ cm^^ (above paragraph), we get F = 
2.73 ± 0.04 and xl = 1-0 {v = 659). When fitting the 
extended component with thermal models, we use the new 
models invoked in XSPEC version 11.0.1 ^ . We can reject 
a CEI model (which describes a plasma in ionization equi- 
librium) with solar abundances. In Fig. 7, we show the 
spectrum of the extended component fitted with a power 
law (left) and APEC model (right). The latter model is 
the new code of thermal plasma including updated atomic 
lines ^ . It is clear that the APEC model over-estimates 
the data in the Fe-K line region. A very low Fe abun- 
dance is required to get an acceptable fit, with the follow- 
ing parameters: 7V//=1.46 x 10^^ cm^^, kT = 4.1 keV, 
Fe = Ni =0 (< 0.16, 2cr), xl = 0.96 {y = 657). While 
the fit is formally acceptable, the derived Nh is incon- 
sistent with that of the SNR. We have subsequently fit- 
ted the extended component with non-equilibrium ioniza- 
tion (NEI) models, which are more appropriate for mod- 
eling SNRs whose age is smaller than the time required to 
reach ionization equilibrium. We use the PSHOCK model 
(Borkowski, Lyerly, and Reynolds 2000), which comprises 
a superposition of components of different ionization ages 
appropriate for a plane-parallel shock. This model is char- 
acterized by the constant electron temperature, Tg, and 
the shock ionization age, Uet (where rig is the postshock 
electron density, and t is the age of the shock). This 
model yields an adequate fit, but with an unusually low 
ionization time-scale and again a very low value of Nh- 
The parameters are: Nh = 1.53 (1.48-1.57) xlO^^ cm"^^ 
kTs = 3.50 (3.28-3.74) keV, n^t = 1.1 (0.8-1.3) xlO^ 
cm^"^ s, and x^ = 0.86 (j^=657). We summarize these 
results in Table 1. When fitting with the NEI model, 
which is characterized by a constant-temperature and a 
single-ionization timescale, we find a similar result to the 
PSHOCK Tuodel. Varying the abundances of Mg, Si, S, Fe, 
and Ni (tied to Fe), the PSHOCK and NEI models yield 
a much higher tiq t (~10"'^^ cm^"^ s), but require unusu- 
ally low abundances. These models mimic the CEI model 
with low abundances. We conclude that while thermal NEI 
models give adequate fits, they either require an unusually 
low Uet (with solar abundances), or alternatively a large 
Uet with unusually low metal abundances. Furthermore, 
the derived Nh from any thermal model is inconsistent 
with that derived for the inner core. To account for a 
change of (5(7V//) > 0.5x10^^ cm~^ ^ J Ue dl, an unrealis- 
tic local electron density of Ue > 10"^ cm~^ is needed for a 
path I ^ 2.4 pc (the extent of the extended component). 
We subsequently examined the spectral variations across 
the SNR. For that purpose, we extracted spectra from the 
brightest knots 2' north of the core, the north-eastern edge 
of the SNR, and an inner circle south-east of the inner core. 
The selected regions are shown in Fig. 8 (North is up. East 
is to the left), and their corresponding spectra (in the 0.5 
- 8.0 keV) are shown in Fig. 9. There is no evidence of 
line emission. All the spectra are well fitted with a power 
law model with Nh ~ 1.9 xlO^^ cm~^ (tied for the three 
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spectra), r=: 2.14, 2.19, and 2.5 for the northern knots, the 
north-eastern edge, and the south-eastern region, respec- 
tively (x^ = 0.50, V = 416). The spectrum of the northern 
quadrant seems to have a harder spectrum that the rest of 
the diffuse emission. This result is in agreement with the 
hardness ratio map shown in Fig. 2. 

4.2. ROSAT PSPC 

The source spectrum of the inner core is extracted from 
a circle of radius R = 40". This defines the extent of the 
plerion seen with ROSAT (Fig. 3). When fitting the PSPC 
spectrum in the 0.1-2.4 keV with a power law model, we 
get Nh = 1.95 (1.35-2.74) x 10^2 cm'^, F = 1.95 (0.5-3.8), 
and xt = 0.95 {ly = 75). The errorbars are large due to the 
poor statistics and the narrow energy band of ROSAT. To 
better constrain the parameters, we fit the ROSAT PSPC 
(0.1-2.4 keV) spectrum simultaneously with the ACIS-S 
(0.5-10 keV) spectrum extracted from the 40"-radius core. 
The power law model yields: Nh = 2.23 (2.22-2.76) xlO^^ 
cm-2, r = 1.85 (1.83-1.87), and xl = 1-1 (t^ = 889). 
The total observed fiux from the inner core is 5xl0~^^ 
ergs cm~2 s~^, which translates to an unabsorbed lumi- 
nosity Lx — 2.8x10^^ erg s^^ (at a distance of 5 kpc). 
The data and fitted model are shown in Fig. 10. 

4.3. Comparison with ASCA archival data 

We here present the analysis of the ASCA archival data 
in order to compare ASCA and Chandra spectra. While 
ASCA has a much poorer angular resolution than Chan- 
dra, its spectral resolution is comparable to that of the 
ACIS-S3 chip. ASCA has the advantage over EXOSAT 
and GINGA in better determining the average photon in- 
dex of the plerion. This was noted by Gallant & Tuffs 
(1998) who found that extrapolating from X-rays to the 
infrared wavelengths using the photon index determined 
with EXOSAT and GINGA yields fluxes lying well below 
the measured ISOCAM values. 

We have extracted the ASCA spectra from a circular 
region of radius i? = 3', centered at the maximum of emis- 
sion (RA (J2000) = 18'' 33™ 33^5, and Dec (J2000) = 
— 10° 34' 3".l). We chose a region of 3' radius since the 
source is 7'. 25 offset from the telescope axis, and a larger 
extraction radius would include events from the edge of the 
detectors. Since G21.5— 0.9 is located in the Galactic bulge 
region, we performed the background subtraction using a 
region from the same field of view. This has the advantage 
of minimizing the contamination of the source spectra by 
the high-energy emission from the Galactic ridge. 

We combine the spectra from the SIS and CIS detec- 
tors, and fit them simultaneously, by introducing a rela- 
tive normalization between the SIS and the CIS spectra. 
A power law model yields an adequate fit with an interstel- 
lar absorption Nh = 2.06 (2.0-2.12) xlO^^ cm^^^ g^ pj^Q. 
ton index F = 1.89 (1.85-1.93), and xl = 0.98 (j/=755). 
The observed flux in the 0.5-10 keV range is 7.6x10"^^ 
ergs cm~2 s~^, which translates to an unabsorbed lumi- 
nosity Lx (0.5-10 keV) = 4.3x10^^ erg s^^ (at a distance 
of 5 kpc). The derived photon index of 1.9 is steeper than 
that derived with EXOSAT (Davelaar et al. 1986), and 
is closer to the value of 2.0 required by Gallant and Tuffs 



(1998) to account for the measured infrared fluxes. The 
small discrepancy might be due to the contamination of 
the ISOCAM fluxes by emission hues or due to subtrac- 
tion of stellar contribution (Y. Gallant, private communi- 
cation). 

5. TIMING 

To search for pulsations in the ASCA data, we extract 
the events from an aperture of radius 3' in the GIS images. 
The GIS data were obtained in the high-bit rate mode, 
with a photon time-of-arrival (TOA) resolution of 0.061 
milliseconds. We perform a power spectral density analy- 
sis on the barycenter corrected photon arrival times using 
the hard band (2-10 keV). No signiflcant pulsations were 
found. The failure to detect pulsations could be partly at- 
tributed to the contamination of the extracted events by 
the plerion. The analysis of the Chandra ACIS-S3 image 
indicates that only ^ 20% of the counts collected in the 
3' ASCA GIS PSF originate from the central 5" radius 
circular region (in the 2-10 keV band). 

We have also searched for pulsations in the Chandra 
HRC data. Slane et al. (2000) analyzed one HRC-I set 
(ObsID I4O6) and obtained an upper limit of 40% on the 
pulsed fraction of the central source. By the time of our 
analysis, there were eleven HRC observations in the pub- 
lic archive. We use only flve of them (Table 2): 2 sets 
acquired with the HRC-S and 3 sets with HRC-I. The rest 
of the HRC observations were observed with large ±Y off- 
sets ^°, that led to considerable spread of the source profile, 
or contained serious processing errors (ObsID 1298). Var- 
ious parameters of the analyzed sets (including the time 
spans Tspan and the effective exposures Toxp) are presented 
in Table 2. All the data were obtained with non-optimal 
focus positions (see parameter SIM_X in Table 2), which 
led to broadening of the core. For each of the 5 sets, a 
short period of time at the start of the exposure was ex- 
cluded from consideration because of an unstable count 
rate. The total time span for all the observations (Sep 
4- Feb 15) amounted to 14,328 ks, with 75.731 ks total 
effective exposure time. 

The HRC-I data were filtered making use of the HRC 
screening algorithm^^. The TOA's were converted to 
barycentric dynamic times (TDB) with the aXhary task 
of CIAO. We extracted the events for the timing analysis 
from the apertures of 15-pixel (« 2") radius centered at 
the core position, the numbers of counts are given in Ta- 
ble 2. The event count rate in set 143 is lower than in 
the other sets because its higher background count rate 
saturated telemetry. 

Because of a wiring error in the HRC detector^^, the 
TOA ascribed to an event in a data file is actually the TOA 
of the following event, which may be a photon event or, 
more likely, a background particle event not telemetered 
to the ground. Therefore, a typical TOA error is equal 
to the average time between the successive events in the 
HRC detector, including the particle events. Making use of 
the Level supplementary data products, which contain 
the total (photon -I- particle) count rates for the whole 
HRC detector, we estimated the mean errors of about 5 

-'^ ''http://asc.harvard.edu/udocs/docs/POG/MPOG/indcx. html 
-'^-'^ http://asc.harvard.edu/cal/Links/Hrc/CIP/filtcr. html 



ms for all the sets except for l^S, for which the error is 
about 1 ms. In principle, the errors could be reduced by 
additional filtering of the extracted events, retaining only 
events such that the time to the next recorded event is 
less than a given value. This filtering, however, strongly 
reduces the number of events available for the analysis, so 
it is not applicable in our case due to the small number of 
photons. Therefore, the timing errors force us to choose a 
minimum period considerably longer than 5 ms (maximum 
frequency smaller than 200 Hz). We chose the frequency 
range 0.01 < / < 50 Hz. 

The presence of several data sets allows us to use the 
"correlated periodicity search" (Pavlov, Zavlin & Triimper 
1999). For each of the data sets, we use the Zf (Rayleigh) 
test (Buccheri et al. 1983). The function Zf{f) has nu- 
merous peaks, with characteristic widths 6f « T'span! o^i^ 
of these peaks may correspond to the pulsar frequency, the 
others are caused by statistical noise. If there is the pulsar 
peak at a frequency /o in one of the data sets, it should also 
be seen in each of the other sets, at frequencies /o + /o AT, 
where /o is the pulsar frequency derivative, and AT is the 
time between the corresponding observations. To identify 
the true pulsar peak (i.e., to find /o and /o), we choose a 
reference data set and examine highest peaks in this set for 
the presence of counterpart peaks in the other sets, for a 
reasonable range of /. If no counterpart peaks are found in 
at least one of the data sets, we have to conclude that the 
peak under examination is due to statistical noise. If coun- 
terpart peaks are found in all the sets, they must satisfy 
some statistical criteria to be considered viable candidates 
for the pulsar peaks (see Koptsevich and Pavlov 2001 for 
a description of the algorithm) . 

For the frequency derivative range, we choose — 1 x 
10^^ < / < Hz s^^, which includes the / values for 
all the known rotation-powered pulsars. All the sets, ex- 
cept for I4O6, meet the requirement |/|maxTspan/2 <C 1, 
which allows us to consider Zi as independent of / for 
each of these sets. To make the computations uniform, we 
divide set I4O6 in two successive subsets, 140(P and 140(f' 
(see Table 2), which satisfy this requirement. 

We applied this approach to the six data sets and found 
no statistically significant periodicity. The series of high- 
est peaks found in all the six sets corresponds to ephemeris 
parameters / = 47.7154362408 Hz and / == -6.60 x 10"^" 
Hz s-i (epoch 51425.25047454 MJD TDB). Combining 
the 6 sets in one, we find a Zf value of 36.6 for these 
/ and /. The probability to find such a peak is close 
to unity, so that the result is statistically insignificant. 
The corresponding pulsed fraction for a sinusoidal signal, 
fp =. {2Zl/Nf/'^ = 16% (where N = 2881 is the total 
number of events), coincides with that measured from the 
extracted light curve (15% ± 5%) and can be considered 
as an upper limit. This is considerably smaller than the 
limit derived by Slane et al. (2000). 

6. DISCUSSION 

We have presented archival X-ray observations of 
G21.5-0.9 with Chandra, ROSAT, and ASCA. The spec- 
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trum of G21.5— 0.9 is non-thermal and best fitted with a 
power law model, whose index steepens by Ar~f from 
the central part of the remnant to the outer compo- 
nent (at ~150"). The total unabsorbed flux from the 
5 arcminute-diameter SNR is Fx (0.5 - 10 kcV) = 1.1 



X 10" 



erg cm 



^ with an 85% contribution from 



the 40" radius core. At a distance of 5 kpc, this trans- 
lates to a total luminosity Lx (0.5 - 10 keV) = 3.3 x 
]^q35 gj-g g-i^ rpj-^g low-surface-brightness extended com- 
ponent discovered with Chandra does not have the limb- 
brightened morphology expected from a SNR blast wave 
or shocked ejecta. It has rather a filled morphology, with 
knots and filaments most evident in the northern quad- 
rant. There is no significant evidence of thermal emission 
lines from this component. Based on the morphology and 
the spectrum of this component, we conclude that it is an 
extension of the inner core, and that the filaments could be 
plerionic 'wisps' as seen in Crab-like plerions. We note that 
XMM-Newton has recently observed G21.5— 0.9 (Warwick 
et al 2000), and the results agree with our findings. In 
this paper, we have investigated in more detail the spec- 
tral variations across the remnant on smaller scales, and 
determined a more stringent upper limit on the pulsed 
fraction. We have also presented for the first time the re- 
sults from the ROSAT and ASCA observations, as a con- 
sistency check for the overall spectral analysis. We show 
that only the non-thermal model fits to the extended com- 
ponent (with Chandra) yield a column density consistent 
with that found for the core, whereas the thermal models 
(CEI or NEI) yield a much lower Nh value. 

In the following, we use our spectral and timing re- 
sults to determine the parameters of the plerion, and a 
'hidden' pulsar. By comparing the Chandra ACIS image 
with the radio image, we conclude that the morphology 
of G21.5— 0.9 could be associated with a pulsar wind and 
possibly jets interacting with the surroundings. The fail- 
ure to detect pulsations could be due to a geometrical ef- 
fect. We estimate the parameters of the putative pulsar in 
G21.5— 0.9 plerion, and compare them with those derived 
for the well-studied Crab and 3C 58 nebulae. Finally, the 
absence of thermal emission in G21.5— 0.9 could be due to 
a low density medium or a cavity in which the supernova 
explosion occurred. 

6.1. Distance 

To estimate the distance to G21.5— 0.9, we regard the 
Nh values found from the ACIS fits, and those from the 
ASCA and ROSAT fits, to describe a reasonable range 
of possible Nh values, (1.8 ~ 2.3) x 10^^ ^m-^. This 
range encompasses the values from all the global power- 
law fits, and excludes the considerably different values re- 
quired for thermal fits. The extinction per unit distance 
in the direction of G21.5— 0.9 can be estimated from the 
contour diagrams given by Lucke (1978): {Eb-v)/D ^ 
0.8 mag kpc"\ Using the relation (Nh/Av) = 1.79 x 
10^^ cm~^ mag~^, which translates into (Nh/Eb-v) = 
5.55 X 10^1 cm-2 mag^^ (Predehl & Schmitt 1995), we 
derive a distance of 4.1-5.2 kpc. This range is consistent 
with the lower estimate on D obtained from 21 cm ab- 
sorption (Becker & Szymkowiak 1981). In the following, 
we assume a distance of 5 kpc and express our results as 
functions of D^ = D/{5 kpc). 



6.2. The putative pulsar 

The morphology and spectrum of the inner core indicate 
the presence of a pulsar powering G21.5— 0.9. The failure 
to detect pulsations could be due to a beaming effect. We 
note that 3C 58 is another Crab-like plerion whose proper- 
ties hint at the presence of a pulsar, yet no pulsar has been 
found there (Torii ct al. 2000). To estimate the spin-down 
energy loss of the pulsar, E, one can use an empirical re- 
lationship connecting E and the X-ray luminosity, Lx, for 
pulsar-powered plerions. For instance, using the formula 
by Seward & Wang (1988), \ogLx = 1.39 logi;- 16.6, 
where Lx is the luminosity of the pulsar plus plerion in 
the (0.2-4) keV band, we get E = 2.7x 10^'^ Dl'^^ erg s^^. 
Using ix = 3 X 10""i;i-22 erg s^^ (Slane et al. 2000) and 
a total luminosity (inner core-|-extended component) of 3.3 
X 10^^ -D5 erg s^^ in the 0.5-10 keV band, we estimate E 

— 5.5 X 10"^^ Dl'^"^ erg s~^. Thus, a plausible estimate for 
the spin-down energy loss is -E37 = _E/(10'^^ erg s^^) ~ 3- 
6. 

We can now use this E and a plausible age of G21.5— 0.9 
to estimate the parameters of the putative pulsar power- 
ing G21.5— 0.9. The age of G21.5— 0.9 is rather uncertain 

— various estimates range from 0.8 to 40 kyr (Salter et al. 
1989, Fiirst et al. 1988, Morsi & Reich 1987). The linear 
size of the X-ray plerion in G21.5— 0.9 is ~ 3.6£'5 pc diam- 
eter, slightly bigger than that in the Crab and 3C 58 nebu- 
lae. The ratio Lx/E for G21.5-0.9 is - 0.006-0.011. This 
is smaller than those for the 1,000-yr-old plerions around 
the Crab pulsar and PSR 0540-69 (0.05), and it is compa- 
rable with Lx/E K, 0.01 for the plerion associated with the 
1,500-yr-old pulsar, PSR 1509-58 in the SNR MSH 15-52 
(Chevalier 2000). On the other hand, the known pulsars 
with £^37 ~ 3-6 have characteristic ages, r, in a range of 
1.5-10 kyr. In the following analysis, we normalize the 
pulsar parameters to T3 = t/(3 kyr). 

Given E and r, the period P of a pulsar can be esti- 
mated as 
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where n is the braking index (P = KP^ 



(1) 



P/[in- 



1)]P] ^ t + To, t is the true age, tq = Po/[in - l)Po], 
Pq and Pq are the initial period and its derivative, and 
/ = 10'*^ /45 g cm^ is the moment of inertia. In this model, 
T Ki t a t ^ To, which is equivalent to P 3> Po- The 
corresponding period derivative is P = P/[{'>^ ^ 1)''"] — 

7.6 X 10-i3[2/(n- l)Y'^lli'^E-:^''^T^^''^. Assuming n = 3 
(magneto-dipolc braking), we can estimate the conven- 
tional magnetic field, P = 3.2 x 10^'^ll^'^R[^{PPf/^ = 

3.2 X lOi^/^/^Pg^/ZVa^^Pj^ff G , where R = lOPio km is 
the neutron star radius. The inferred magnetic field is 
comparable to those of very young radio pulsars. 

The X-ray luminosities, L^ of the radio pulsars de- 
tected with ROSAT (0.1-2.4 keV band) show strong cor- 
relation with the value of E, which allows us to es- 
timate the X-ray luminosity of the putative pulsar in 
G21.5— 0.9. Using the relationship derived by Ogelman 



(1995), L 



(p) 

X 



6.6 X W^{Bi2/P'^f-'^ erg s-\ which is 



equivalent to L^ 



1.5 X 10 i6£;i.35 gj.g g 1^ ^Q obtain 
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1.3 X IQ^'^El^^ erg s ^. The formula suggested 



by Becker & Triimper (1997), L 

(p) 

X 



(p) _ 



X 



10 E, gives ap- 



1 X lO^^i;,, erg s-i 



proximately the same value, L 
We can now compare these values with an observed upper 
limit on the pulsar flux. To estimate a conservative upper 
limit, we measure the flux in a 0".5-radius aperture placed 
at the brightest part of the core (see Fig. 5). We normal- 
ize the PSF to match the number of counts in the central 
0".5 bin. We then estimate the ratio of counts within 0".5 
to counts within 1" radius and scale their corresponding 
fluxes accordingly. Using the estimated flux of3.5xl0^^^ 
ergs cm^^ s^^ from the 1" central core (§ 4.1), we estimate 
an unabsorbed flux of 2.4xl0^^^ergs cm^^ s^^ from the 
point source. This translates to an upper limit on the flux 
from a hidden pulsar of 1.1 x 10^^^ ergs cm~^ s^^ in the 
0.1-2.4 keV range (assuming F — 1.5). If the pulsar's X- 
ray emission were isotropic, the corresponding upper limit 

on the pulsar luminosity would be U~^' < 3.5 x 10'^^ Df 
erg s~^, an order of magnitude lower than the above pre- 
dictions based on the E estimates. This indicates that the 
X-ray radiation of the pulsar is strongly beamed, and the 
orientation of the pulsar's rotation and magnetic axes is 
such that the rotating beam does not come close to the 
line of sight, in agreement with the lack of observable pul- 
sations. 

6.3. Pulsar Wind Model 

In the Kennel and Coroniti (1984; hereafter KC 84) 
model, a pulsar injects high-energy particles in the form 
of a relativistic particle wind. Conflnement of the pulsar 
wind by its surroundings causes the outflow to deceler- 
ate at a shock, forming a synchrotron nebula. The radius 
of the shock, i?s, is given by equating the pressure of the 
pulsar's wind, E/{AitR'^c), with the pressure in the nebula, 
P„. Beyond this radius, a non- relativistic flow transports 
the plasma from the shock region to the edge of the neb- 
ula. The size of the nebula, i?„ is related to the shock 
radius, Rs as: Rn/Rs ^ Vv^! where a is the so-called 
pulsar wind magnetization parameter deflned as the ratio 
of the Poynting flux to the particles flux (KC 84). Self 
consistent models for the Crab require that the wind is 
particle dominated, i.e. a is small (< 10~^; KC 84). From 
the ROSAT PSPC and Chandra ACIS-S spectra, we esti- 
mate an equipartition magnetic field, B^q = 1.8xl0~ G 
in the 40" radius core. This yields a pressure in the in- 
ner 40" radius core of B'^/Stt ^ 1.29x10^^ erg cm""^, and 
therefore Rs — 0.046 £^37^ pc. Using the estimated S37 ^ 
3-6 (above subsection), we get Rs ^ 0.08-0.11 pc; which 
corresponds to an angular size of ^ 3"-5" (at 5 kpc). Us- 
ing the nebular size i?„ ^ 150", we estimate the pulsar 
wind parameter, a ~ {Rg/Rn)'^ ~ (4-11) x 10~^, indi- 
cating a particle dominated wind — like the Crab and the 
other Crab-like nebulae. In Table 3, we summarize the 
inferred parameters of the putative pulsar in G21.5— 0.9, 
in comparison with the Crab pulsar and the putative pul- 
sar in 3C 58. When fitting the spectra in circular rings of 
G21.5— 0.9, we have found that the photon index steepens 
away from the center all the way out to a radius of 150" — a 
result expected from synchrotron losses. Chevalier (2000) 
has developed a one-zone, shocked wind model, to model 



the X-ray luminosity of pulsar nebulae. This model differs 
from the model of KC 84 in that the synchrotron burn- 
off leads to a steepening of the frequency spectral break 
by 0.5, whereas it could be greater than 0.5 in the KC 84 
model. For G21.5— 0.9, the photon index, F, varies from 
1.5 (in the innermost 5" core) to ^ 2.7 (in the outer diffuse 
component); i.e. the net change of ^ 1.2 is inconsistent 
with the Chevalier model. 

6.4. Correlation with the radio image 

The high-resolution 22.3 GHz radio image of G21.5— 0.9, 
obtained with the NRG Millimeter Array, shows a diffuse 
centrally peaked component plus axisymmetric filaments 
(Fiirst et al. 1988). The authors suggest that the sym- 
metry in the lobes is evidence of precessing jets from a 
hidden pulsar. In order to compare the radio image with 
our X-ray image, we scale the ACIS-S image to that of 
the radio image (Fig. 7 in Fiirst et al. 1988) and show the 
resulting image in Fig. 11. The left panel shows the 40" 
radius inner core of G21.5— 0.9, binned to the same scale, 
smoothed with a Gaussian with cr = 3" , and with the con- 
tours overlayed. The right panel shows the radio image 
of the core with the jets model proposed by Fiirst et al. 
(1988). It is clear that the X-ray trough in the north-west 
is strongly correlated with that of the radio image. Our 
result therefore supports the proposed picture. It is inter- 
esting that the shock radius {Rs ^ 0.1 pc; § 6.3) coincides 
with the axisymmetric filaments seen in the radio (Fiirst et 
al. 1988). Therefore, those structures could be due to the 
deceleration of the pulsar wind (or jets) confined by the 
nebular pressure. They also found that the polarization 
varied from the center to the outer boundary, from almost 
zero to 20%-30%. They conclude that this is explained by 
a change in the magnetic field geometry, with little contri- 
bution from a radial field from the inner regions because 
of geometrical effects. 



6.5. Whe 



the SNR shell? 



To date, we know ^ 10 Crab-like SNRs which lack a 
SNR shell, which would be expected from the shocked am- 
bient matter or shocked ejecta. The lack of SNR shells is 
still a puzzling issue, and many deep searches for shells 
around the Crab and 3C 58 have failed to find them. It 
has been suggested that the absence of a shell could be 
due to a propagation in a low density medium (in which 
case the shell would be very faint), or to a low efficiency 
in accelerating particles and amplifying the magnetic field 
(Reynolds & AUer 1985). A deep search in the radio for 
G21.5— 0.9 shell is needed, since all current radio obser- 
vations have been performed to study mainly the bright 
inner core. Slane et al. (2000) suggested that the extended 
component is associated with the blast wave of the SN ex- 
plosion, or from the shocked ejecta. Here we rule out ther- 
mal CEI models. Fitting with NEI models, which are more 
realistic to model thermal emission from young SNRs, in- 
dicates unusually low abundances or unusually low ioniza- 
tion time-scales. Furthermore, the derived Nh is inconsis- 
tent with that derived for the SNR. We find no evidence 
of line emission in any part of the remnant, leading us to 
conclude that the spectrum is dominated by non-thermal 
emission. 

A growing number of SNRs show evidence of non- 
thermal emission from their shells (e.g. SN 1006). The 



emission mechanism has been attributed to synchrotron 
radiation from high-energy particles accelerated at the 
SNR shock (Allen, Gotthelf, & Petre 1999). Reynolds and 
Chevalier (1981) have noted that the non-thermal spectra 
should not be described by a power law, but rather a power 
law with a roll-off at high energies. If the non-thermal 
emission in the extended component of G21.5— 0.9 is in- 
deed due to synchrotron radiation from highly energetic 
particles accelerated at the SNR shock, then one would 
expect it to be fitted with the cut-off or escape models 
recently included in XSPEC vll.O, which are more appro- 
priate for modeling non-thermal emission from SNR shells 
(Reynolds & Keohane 1999, Dyer et al. 2000). These mod- 
els require knowledge of the radio flux density and the 
radio spectral index. With the absence of a radio coun- 
terpart, we fit the X-ray spectrum with the sresc model 
leaving these parameters free. This model does not pro- 
vide a good fit {xi = 1.36, v — 670). This again argues 
against a shell interpretation for the extended component. 
The non-thermal emission in the extended component 
is most likely due to synchrotron emission from particles 
accelerated and injected by the 'hidden' pulsar. Using the 
power law model parameters of the extended component 
(Table 1, 1st row), and assuming equipartition of particles 
and magnetic field, we estimate an equipartition field of 
1.8x 10""* Gauss. The lifetime of the synchrotron emitting 
electrons is only 60 E~1fy BZ\'^ years; where E^^keV and 
i?_4 are the photon's energy in units of 1 keV and the 
magnetic field in units of 10""* G, respectively. The syn- 
chrotron lifetime is much smaller than the estimated age 
of the SNR — a result also found in the Crab and Crab- like 
nebulae. The electrons can propagate out to a radius of 
3 pc if they are traveling at a speed of 0.15 c E^^-^^y Bif. 
The enhancement of brightness in the northwest quadrant 
could be due to the propagation into a denser medium 
(as evident in the public IRAS images of G21.5— 0.9). To 
get an upper limit on the density of the ambient medium, 
we fit the extended component with a power law plus a 
thermal model (APEC). We freeze the power law model 
parameters to their best value (Table 1), and determine an 
upper limit on the emission measure (EM) of the thermal 
component. We get EM < l.lxlQ-* (at the 90% con- 
fidence level). Using EM = 10"^'' {4ttD'^)-'^ JneUHdV, 
where rig and uh are the postshock electron and Hydrogen 
density, respectively, and V is the volume of the extended 

component, we estimate {nf-nn) < 0.07 D^ cm^''. 
Assuming a strong shock with compression ratio 4, the 

ambient density, no, is < 0.02 Dc^ cm~'^. 

G21.5— 0.9 is perhaps most similar to the plerionic SNR 
3C 58, although there are some differences. Both plerions 
have a morphology hinting at the presence of a pulsar, 
yet undetected. They both have a low-frequency spectral 
break (at 50 GHz), and thus belong to the class of ple- 
rions of the 'second kind' (Salvati et al. 1998). Woltjer 
et al. (1997) suggest that these plerions are powered by 
short-lived pulsars with either very small braking indices 
or pulsars that would undergo a phase change in the pul- 
sar's energy output. Both are young and exhibit a bright 
inner core plus a fainter extended component whose spec- 
tra are dominated by non-thermal emission. The extended 
component in 3C 58 shows however a small flux (6% of the 

^^http://cac. psu.edu/bcatnic/Cluster/Lionx/lionx. html 



total) of thermal nature, which was attributed to the ex- 
pansion of the outer rim into inner ejecta core (Bocchino et 
al. 2001). G21.5— 0.9 is also similar to the Crab nebula in 
that it has a centrally-filled morphology and a non-thermal 
spectrum. The softening of the spectral index away from 
the center in G21.5— 0.9 is similar to that in the Crab 
nebula (Willingale et al. 2001). In spite of these simi- 
larities, G21.5— 0.9 remains a unique plerion. Its X-ray 
brightness profile (core-|-plateau) and its morphology are 
puzzling, and can not be explained by diffusion models 
such as Becker (1992). The detection of axisymmetric fil- 
aments in the radio and X-rays hint at a complex magnetic 
field geometry and/or jet structures (Warwick et al. 2000). 
The most unusual property of G21.5— 0.9 plerion lies 
in the radio regime. Plerions are usually bigger in size 
at lower energies, since the high-energy photons have a 
shorter lifetime than the low-energy photons. There- 
fore, one would expect to see a bigger plerion in the ra- 
dio. This is the case for the Crab nebula. In 3C 58, 
there is a close correspondance between the radio mor- 
phology of the outer nebula and the X-ray map at soft 
energies (Bocchino et al. 2001). For G21.5— 0.9, only 
the inner 40" radius core has been detected in the ra- 
dio. Slane et al. (2000) have put an upper limit on the 
1 GHz surface brightness of the extended component of 
E = 4x 10-21 w m-2 Hz-i sr"!. Gaensler (2000, private 
communication) has recently found a radio knot coincident 
with the hard X-ray brightest knot/filament found in the 
northern region. A deep observation of G21.5— 0.9 in the 
radio and at the short wavelengths (currently conducted 
with the VLA) will determine the spectral characteristics 
of this knot and will allow us to search for the radio coun- 
terpart of the extended component and the missing SNR 
shell. 
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Fig. 1.— The Chandra ACIS-S images of G21. 5-0.9 at 0.1-6.0 keV (top) and 6.0-10.0 keV (bottom). The total exposure time is 71.6 ksec. 
The images are smoothed with a Gaussian with ct = 3". Coordinates are J2000. 
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Fig. 2. — Hardness ratio map of G21.5— 0.9 observed with Chandra ACIS-S. The ratio is 2.4-10 keV over 0.5-2.4 keV. The image is smoothed 
with a Gaussian with cr = 3". 
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Fig. 3.— The ROSAT HRI contours of G21.5-0.9 showing the plerion {_R~40") overlayed on the soft band image obtained with ACIS-S3. 
The point-Uke source southwest of the plerion coincides with an emission-line star, SS 397. Both the HRI and ACIS-S3 images are smoothed 
with a Gaussian with u = 6" and contours are displayed in a log scale. 
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Fig. 4.— The HRC-I image of the inner core of G21. 5-0.9. See S3.3 for details. 
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Fig. 5. — The radial profile of G21.5— 0.9 using HRC and ACIS. The HRC profile extends out only to a radius ~ 40" since the background 
dominates at larger radii. The inset shows the HRC profile with 0".25 bins and the expected point-spread function (PSF, dashed line). 
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160 



Fig. 6. — The Power Law photon index F as a function of radius (in arcseconds). Zero radius corresponds to the brightest emission region 
of the plerion, centered at RA = 18'^33™33''.50, Dec = — 10°34'6".6 (J2000). We show the power law index of the inner 40" radius core (in 
radial increments of 5") and the extended component (50"-150"), for Nh = 2.24 X 10^^ cm~^. 
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Fig. 7. — The spectra of the extended component (50"-150") fitted with a power law model (left) and the APEC thermal model assuming 
solar abundances (right). We show that while a non-thermal model fits this component adequately, the thermal model fails. The parameters 
of the best fit (power law) are summarized in Table 1. 
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Fig. 8. — The Chandra ACIS-S (0.5-10 keV) image of G21.5-0.9 showing the inner core (i?~40") and the diffuse component extending 
out to a radius of _R~2'.5. We also show (in green) the regions selected for the spectral analysis: the north-eastern knots, the brightest 
north-western edge, and the south-eastern region. 
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Fig. 9. — The data and the fitted power law models for selected regions from the extended component, shown in Fig. 8. The blue, 
and red colors correspond to the bright north-eastern knots, the north-western edge, and the south-eastern region, respectively. 



green. 
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Fig. 10.— The PSPC and ACIS-S3 spectra of the inner 40" core fitted with a power law model (§4.2). 
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Fig. 11. — The ACIS-S and 22.3 GHz images of G21.5-0.9, binned to the same scale (B1950 coordinates). The morphology of the plerion 
in X-rays follows the enhancements seen in the filamentary radio structures. 
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Table 1 
Spectral parameters of the fits to the low-surface-brightness extended component using ACIS-S3 

OBSERVATIONS. THE RANGES ARE AT THE 90% CONFIDENCE LEVEL. 



Model 


iVH^(xlO^^ cni-^) 


Model Parameter 


Norm. (xlO-^) 


xii^) 


Power Law 


1.83 (1.77-1.90) 


r = 2.36 (2.30-2.43) 


4.36 (4.04-4.78) 


1.0 (658) 


APEC'^ 


1.63 (1.58-1.69) 


kT = 3.33 (3.14-3.54) 


8.22 (7.89-8.59) 


1.15 


P SHOCK 


1.53 (1.48-1.57) 


fcT = 3.50 (3.28-3.74) keV 
t'== 1.1 (0.8-1. 3)xl09 cm-3 s 


9.72 (9.23-10.2) 


0.86 (657) 



"Best fit values are indicated. When freezing Njj to 2.2 
thermal models are rejected (x2"^l-5) 



X 10^ 



^ the power law model gives T = 2.73 ± 0.04 (xl = 0.98; u = 659), and 



''Assuming solar abundances (ittp://hea-www. harvard.edu/APEC) 

'^The ionization timescale, rtet; where ne is the postshock electron density, and t is the age of the shock 
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Table 2 
Chandra/URC observational sets 



ObsID^ 


Detector 


Date 


Epoch'^ 


rp C 

-'- span 


rp d 

-'- exp 


SIMx'' 


N' 








MJD 


S 


S 


mm 




143 


HRC-I 


Sep 4, 1999 


51425.339 


15 243 


15 243 


-1.0160 


179 


1242 


HRC-I 


Sep 4, 1999 


51425.524 


14 855 


14 855 


-1.0160 


555 


147 


HRC-S 


Sep 5, 1999 


51426.625 


9 747 


9 747 


-1.3705 


496 


I4O6'' 


HRC-I 


Oct 25, 1999 




27 572 


25 621 


-1.0389 




146 


HRC-S 


Feb 16, 2000 


51591.026 


10 334 


10 265 


-1.4282 


518 


UOff" 






51476.316 


13 797 






543 


140^ 






51476.475 


13 726 






590 



''the ObsID I4O6 set was divided in two parts, see explanations in Section 5 

''epoch refers to the central time of the exposure 

'^time span 

■^effective exposure time 

■^focus position along the telescope axis; the optimal positions arc —1.04 mm for HRC-I, and —1.53 mm for HRC-S 

•'^ number of counts used in the analysis 
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Table 3 
Comparison of the parameters of the putative pulsar in G21.5— 0.9 with the Crab pulsar and the putative 

PULSAR IN 3C 58 (Torii et al. 2000). 





G21. 5-0.9 


Crab 


3C58 


Distance (kpc) 


5 


2 


3.2 


Lx^(1035 erg s'l) 


3.3Z)i 


210 


0.2 


E (1037 erg s-i) 


3-6 


47 


0.4 


ab 


(4-ll)xlO~4 


3x10-3 


- 6x10-3 


P (ms)'= 


144r3-°-54-,°-^ 


33 


440 


Bo (1013 Gauss) 


1.1 T^^E^r^-"" 


0.4 


6 



"X-ray luminosity of the plcrion in the 0.5-10 keV range 

''Pulsar wind parameter defined as the ratio of the Poynting flux to the kinetic energy flux (Kennel and Coroniti 1984) 

"^ The following notations are used: rs is the age of the SNR in units of 3 kyr; B37 is the spin-down energy loss in units of lO'^'^ erg s^^, D5 is 
the distance in units of 5 kpc 



